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ABSTRACT 

The effect of 40 keV electron irradiation on a-Si:H p-i-n single-junction solar cells was investigated using 
measured and simulated dark J-V characteristics. EPRI-AMPS and PC-1 D simulators were explored for use in 
the studies. The EPRI-AMPS simulator was employed and simulator parameters selected to produce agreement 
with measured J-V characteristics. Three current mechanisms were evident in the measured dark J-V 
characteristics after electron irradiation, namely, injection, shunting and a term of the form CV m . Using a single 
discrete defect state level at the center of the band gap, good agreement was achieved between measured and 
simulated J-V characteristics in the forward-bias voltage region where the dark current density was dominated by 
injection. The current mechanism of the form CV m was removed by annealing for two hours at 140 °C. 
Subsequent irradiation restored the CV m current mechanism and it was removed by a second anneal. Some 
evidence of the CV m term is present in device simulations with a higher level of discrete density of states located 
at the center of the bandgap. 


INTRODUCTION 

Thin-film a-Si:H solar cells are prime candidates for space applications because investigations show that they 
have good radiation resistance [1], The fabrication technology may be automated and has the potential for solar 
arrays with low cost per watt. The a-Si:H-based solar cell technology permits deposition on ultra-thin flexible 
substrates, thus offering high specific power density. To take advantage of the high specific power density, cells 
must be flown without a cover glass, thereby requiring that the nude device structure have good radiation 
resistance. In order to predict the radiation resistance of solar cells in space environments, particle irradiations 
are generally carried out in the laboratory with monoenergetic particle beams at normal incidence. In space 
environments there is a wide energy distribution of particles with a range of incident angles. Hence, the 
fundamental radiation resistance mechanisms must be understood if laboratory measurements are to be used to 
predict solar cell performance in a variety of space environments. 

The goal of our program is to acquire a detailed understanding of the effect of particle irradiation on thin-film solar 
cells. We have carried out studies on the effect of MeV proton irradiation on single and multi-junction a-Si:H cells; 
dark current mechanisms have been identified in both virgin and irradiated cells [2], The proton irradiation studies 
have shown that dark J-V characteristics are orders of magnitude more sensitive to radiation effects than light J-V 
characteristics. Woodyard and Landis reviewed the published literature on electron irradiation studies on thin- 
film devices [3]; materials and device structures have improved since the review was published. Cells are thinner 
and more efficient than those studied in earlier work. It is important to study and characterize the effect of particle 
irradiation on state-of-the-art cells if current technology is to be advanced. The AE8MAX model [4] predicts a ten- 
year integrated electron flux of 1.7E16 cm' 2 for electrons in the 50 keV to lOMeV range incident on a surface in 
geosynchronous orbit. The model does not include flux from solar flares and fluxes are questionable at the lower 
electron energies. The flux of electrons in the model increases at lower energies. It is for this reason that our 
initial electron irradiation work was carried out with electron beam energies of 40keV and fluences of about 1 E 1 7 
cm' 2 . 


1 This work was supported under NASA Grant NAG 3-2180. 
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A single-junction a-Si:H p-i-n device structure was selected for this work. Single-junction p-i-n cells without back 
reflectors, as opposed to multi-junction cells with back reflectors, were studied for two reasons. The single- 
junction p-i-n device structure is the basic structure employed in multi-junction a-Si:H solar cells. An 
understanding of the single-junction a-Si:H p-i-n device structure may make it possible to predict the performance 
of multi-junction device structures. Additionally, the analysis of measured dark J-V characteristics of a single- 
junction p-i-n device structure is considerably easier than those of multi-junction devices structures. 

The study of the fundamental radiation resistance mechanisms in devices is complex. It is for this reason that 
both J-V measurements of irradiated cells and device simulations of J-V characteristics were employed in this 
work. In order to gain experience with available device simulators in calculating dark J-V characteristics, two 
simulators were used to study a two-sided step p-n junction. The free electron densities, electric field and dark J- 
V characteristics of a two-sided step p-n junction can be calculated using an analytical model [5], The simulators 
selected were EPRI-AMPS [6] and PC-1D [7], Parameters were selected for each simulator to produce identical 
electron densities, electric field and dark J-V characteristics. The premise for this approach was, if two simulators 
developed independently and employing different algorithms could produce the same p-n junction characteristics, 
and the characteristics agree with an analytical model, then simulations of p-i-n structures, which cannot be 
modeled analytically, may have some validity. Agreement between the simulators and analytical would also 
insure that we had an adequate understanding of the simulators to use them effectively. 

Following irradiations with 40 keV electrons, dark J-V characteristics were measured. Simulations of dark J-V 
characteristics of p-i-n solar cells were carried out and parameters selected to fit the simulations to the 
measurements. The parameters that resulted in good fits to the measured J-V characteristics in the various 
forward-bias ranges were used to develop an understanding of current mechanisms in both virgin and irradiated 
a-Si:H p-i-n solar cells. 


p-n JUNCTION SIMULATIONS 

A p-n junction device structure with Si material parameters was investigated with both PC-1D and EPRI-AMPS 
simulators and the results compared with each other and an analytical model for a two-sided p-n step junction. 

The p-type and n-type layers were each 5000 pm thick. A doping level of 1 El 8 cm’ 3 was used for both acceptors 
in the p-layer and donors in the n-layer. For simplicity, a mobility of 1000 cmW 1 was chosen for both electrons 
and holes. In working with the simulators, we learned there are two major differences between the simulators. 
First, PC-1D is a transient simulator, and EPRI-AMPS is a steady-state simulator. Therefore, stepping times for 
sweeping device voltages must be considered when using the PC-1D simulator. Secondly, the PC-1D simulator 
uses a lifetime model and the EPRI-AMPS simulator employs an electronic density of states, DOS, model. In 
order to calculate free electron carrier densities from lifetimes with the EPRI-AMPS simulator, a lifetime 
subroutine is provided by the developer for calculating a DOS from carrier lifetimes. The PC-1 D and EPRI-AMPS 
simulators produced vastly different free electron densities with the lifetime approach. It was concluded that the 
EPRI-AMPS simulator does not produce correct free electron densities when the lifetime approach is employed. 
The next approach was to use the DOS feature of EPRI-AMPS and to calculate a lifetime using a model that 
relates the DOS to a lifetime that could be used in the PC-1D simulator. 

Figure 1 shows the simulated and calculated free electron concentration in the p-n junction under a forward bias 
of 0.5 V. The front of the device is located at a positon of 0 pm and the metallurgical junction at 5,000 pm. The 
front p-type layer is within positions of 0 to 5,000 pm and the rear n-type layer is within positions of 5,000 to 
10,000 pm. The front and rear contacts in the simulators are located at 0 and 10,000 pm, respectively. Electrons 
are minority carriers in the p-type layer and the free carrier densities are about 1 E2 cm" 3 near the front contact of 
the device. Electrons are majority carriers in the n-type layer and have free carrier densities of about 1 El 8 cm 3 
near the back contact. Figure 1 shows that both PC-1D and EPRI-AMPS simulators produce a diffusion electron 
current density due to electron injection at the metallurgical junction. Electron injection is evidenced by the 
exponential decrease in the free electron carrier densities from the metallurgical junction forward into the p-type 
layer. The PC-1D simulator free electron density versus position was calculated using a Shockley Read Flail, 

SRH, recombination lifetime of 10 ps. The band-to-band recombination coefficient was set to an unrealistically 
low number for simplicity so that the recombination is entirely dominated by SRH recombination. The EPRI- 
AMPS simulator free electron density versus position shown in Figure 1 was calculated with a DOS distribution. 
The DOS distribution used was constant both with respect to position in the device as well as with respect to 
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Figure 1 . Simulated free electron density of a p-n junction with a forward bias of 0.5 V. 


energy in sub-band-gap region of the energy band. The DOS used in the EPRI-AMPS simulation, n E , was 1 E14 
cm^eV' 1 ; the cross section for the states, o t , was IE-16 cm 2 . The thickness of the layers, 5,000 pm, was chosen 
to insure that all recombination occurs in the bulk of the material and not at the contacts. 

The exponential decay of minority carriers in the p-type layer can be fit by an exponential with two parameters, 
namely, the injected excess electron minority carrier density, An p =1E10 cm’ 3 at 5,000 pm, and an electron 
diffusion length, L n =1 60 pm. The fit of the excess electron carrier density, 8n(x p ), to the exponential behavior in 
the p-type layer with x p <5000 pm is shown as the solid line in Figure 1 ; it was obtained using: 

8n(x p )=An p exp[(x p -5000 pm)/L n ] . (1) 

The electron diffusion length can be calculated from the electron mobility, p n , and lifetime, x n , using: 

L n =sqrt [(kT/q)p n x n ) (2) 

where k is the Boltzman constant, T the temperature, q the electronic charge, and x n the electron recombination 
lifetime. Substituting p n =1000 cmW 1 and x n =10 ps into Equation (2) yields L n = 1 60 pm. Hence the lifetime 
model employed in the PC-1 D simulation is in agreement with the analytical theory. 

An equivalent DOS lifetime for comparison with the PC-1 D simulator was calculated from the parameters used in 
the EPRI-AMPS DOS simulation using [8]: 


x n =( n E AEo t Vt)' 1 (3) 

where AE is the energy width of the DOS distribution for states serving as recombination centers and v t is the 
carrier thermal velocity. Using n E =1 El 4 cm' 3 eV _1 , AE=0.1 eV, o t =1E-16 cm 2 and v t =1E7 cm/s, yields a lifetime of 
100 ps, which is ten times greater than the lifetime used in the PC-1 D simulation. The free electron density 
shown in Figure 1 is in agreement with the PC-1 D simulation, and the two simulations are in agreement with 
analytical theory for a two-sided step p-n junction. However, the discrepancy in the lifetimes must be resolved if 
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Figure 2. Simulated electric field of p-n junction near the metallurgical junction with a 
forward bias of 0.5 V. 


one is to have confidence in applying the simulators to the analysis of devices. We will explore the discrepancy in 
future work. 

The position dependence of the electric field was simulated with the two simulators. The results are shown in 
Figure 2. The electric field has the largest value at the metallurgical junction and the depletion width is about 0.05 
pm for both simulations. The step like behavior in the results for the PC-1 D simulation is due to rounding-off 
errors. The electric fields predicted by both simulators are in agreement with the analytical model. 

The simulated dark J-V characteristics of the p-n junction for the EPRI-AMPS and PC-1 D simulators are shown in 
Figure 3. Also shown is the dark J-V characteristic calculated with the analytical model: 

J(V)=J 0 [exp(qV/kT)-1] where J 0 =(kT)( p p p n /L p + p n n p /L n )j. (4) 

Substituting the values for the Boltzman constant, temperature, mobilities, minority carrier densities and diffusion 
lengths and into Equation (4) results in the straight line shown in Figure 3. There is good agreement between the 
analytical model and the simulated results in the 0.2 to 0.8 V range of forward biases. As mentioned above, the 
PC-1 D simulator is a transient simulator, and for this reason, current densities at low voltages with depend on the 
device capacitance and voltage step rate. The effect of varying the step rate of the forward-bias voltage on the 
dark J-V characteristics is shown in Figure 4. The current charges the junction capacitance at low biases. For 
these simulations, the p-n junction capacitance is 0.2 pF/cm 2 . When the voltage is ramped at 10 mV/s, which 
corresponds to 1 s/step, the current at voltages less than 200 mV is 2e-9 A/cm 2 . Reducing the step rate to 100 
s/step results in the forward-biased current decreasing by an order of magnitude. A 1000 s/step decreases the 
current further. The dark J-V characteristic approaches the characteristic that would be expected for the 
analytical model. 
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SINGLE-JUNCTION p-i-n SOLAR CELL MEASUREMENTS AND SIMULATIONS 

Single junction p-i-n a-Si:H solar cells were deposited on bare stainless steel substrates without the use of 
textured back reflectors. An n-type layer about 15 nm thick was deposited on the bare stainless steel substrate 
followed by the deposition of a high quality a-Si:H intrinsic layer about 200 nm thick. An 8 nm thick p-type layer 
was then deposited on the i-type layer. All a-Si:H layers were deposited using radio-frequency plasma enhanced 



Voltage (V) 

Figure 3. EPRI-AMPS and PC-1 D simulated dark J-V characteristic for a two-sided-step 
p-n junction. 



Voltage (V) 

Figure 4. PC-1 D simulated dark J-V characteristics for a two-sided step p-n junction with 
three voltage step rates. 
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chemical vapor deposition. Indium tin oxide, ITO, was evaporated through a mask with an area of 0.25 cm 2 . The 
ITO serves as the top contact and as a % wavelength anti-reflective coating. Finally, Au grids were evaporated 
through a mask on top of the ITO. Since the cells do not employ a light scattering back reflector, the initial 
efficiency of the devices is about 6%. Initial dark J-V characteristics were measured after the virgin cell was 
annealed for two hours at 200 °C in a vacuum system at about IE-6 Torr. A Cambridge S250 MK3 scanning 
electron microscope was used to irradiate the single-junction p-i-n solar cells with 40 keV electrons. The electron 
current was about 0.6 pA and the beam scanned over an area of about 0.55 x 0.72 cm 2 . Irradiation times of five 
hours were necessary to obtain a fluence of about 1 El 7 cm' 2 . Following the irradiation, dark J-V characteristics 
were measured again. The cells were annealed a second time. A subsequent irradiation and annealing cycle 
was conducted. Dark J-V characteristics were measured between each irradiation and annealing step. The dark 
J-V characteristics were measured with a computer-interfaced Keithley Instruments model 236 source 
measurement unit that is programmed to source the bias voltage in 10 mV increments and measure the current at 
1 1 1 separate voltages in the range of 0 to 1.10V. The instrument and software are configured to make 320 
current measurements at each source-bias voltage and calculate the average and standard deviation of the 
current. It took three to four minutes to measure a single J-V characteristic. The standard deviation for each of 
the measurements was calculated and is greater than 100 for the measurements reported in this work. This 
insures that the average current is representative of the device and not a transient current due to junction 
capacitance. 

The energy deposition of the irradiating 40 keV electrons is assumed to be uniform in the cells studied in this 
work. While it is possible to use particle-solid interaction models to determine the spatial dependence of the 
electron-hole generation rate produced by the irradiation, our interest is in determining the role of the material 
properties on the dark J-V characteristics. Since the energy deposition is assumed to be uniform, simulated 
material properties were chosen to be uniform in the intrinsic layer with respect to position. This assumes that 
radiation induced defects do not migrate from the position were they are created. 

The measured dark J-V characteristics of the p-i-n solar cell before irradiation are shown by the dots in Figure 5. 
The dark J-V characteristic measured after annealing at 200 °C for 2 hours is labeled all. A second J-V 
characteristic, al 2, was measured immediately following the all. The currents for both the al 1 and the a12 
measurements are linear on the semi-log graph over a range of 1 E-1 0 to 1 E-4 A/cm 2 . The data show that the 
measurement of al 1 influenced the measurement of a12. Subsequent measurements of the dark J-V 
characteristics are similar to the a12 measurements, thereby showing saturation of the effect that produced the 
a12 dark J-V characteristic. Flence, the cell was annealed before each irradiation. EPRI-AMPS simulations were 
carried out with various material parameters in an effort to obtain a good fit between the measured simulated dark 
J-V characteristics. The most influential parameter that affected the dark J-V characteristics was found to be the 
DOS near the center of the band gap. A good fit was obtained with a DOS made up of a discrete energy level 
with an energywidth of 0.1 eV and located at the center of the band gap. Equal numbers of acceptor and donor 
like states were used for the discrete energy states. The results of simulations with various densities of discrete 
energy states are plotted in Figure 5 as solid lines. The density of discrete energy states was increased from 0 to 
1 El 7 cm' 3 eV' 1 . The simulated results can be fitted with an expression of the form: 

J(V)=J 0 [exp(qV/nkT)-1] (5) 

where n is the material quality factor. As the density of discrete energy states at the center of the band gap is 
increased, J 0 and n increase. The curve labeled 0 corresponds to a zero discrete energy state density at the 
center of the sub-band-gap. The recombination for a simulation with no discrete energy states is dominated by 
recombination from the conduction band-tail states to the valence band-tail states. Figure 5 shows that the J 0 and 
n terms in Equation (5) for the curve labeled 0 are significantly less than for the measured all data. The curve 
labeled 5E15 represents a simulation with a discrete energy state density of 5E15 cm' 3 eV" 1 at the center of the 
band gap. The J 0 and n terms for this simulation are in agreement with the al 1 data over 6 orders of magnitude. 
The curve labeled 1 E 1 6 shows the results for a simulation with a discrete energy state density of 1 E16 cm' 3 eV' 1 
located at the center of the band gap. This simulation is in agreement with the a12 data over 6 orders of 
magnitude as well. The deviation in linearity between the simulated and measured dark J-V characteristics above 
0.8 V is believed to be due to series resistance in the measured cells. The series resistance results from factors 
external to the p-i-n structure, such as contact resistance and resistance in the ITO layer. Future simulations 
including these factors will be carried out to obtain a better fit of the dark J-V characteristics above 0.8 V. 
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Figure 5. Measured and EPRI-AMPS simulated dark J-V characteristics for an a-Si:H 
p-i-n solar cell before irradiation. 


Figure 5 suggests that the difference between the al 1 and a12 J-V characteristics is due to an increase in the 
sub-band-gap energy state density from 5E15 cm' 3 eV" 1 to 1 E16 cm" 3 eV' 1 . 

A dark J-V characteristic of a p-i-n solar cell after irradiation with 40 keV electrons to a fluence of 1 El 7 cm" 2 is 
shown in Figure 6 by the curve made up of dots and labeled rl . For comparison purposes, the measured dark J- 
V characteristic before irradiation is also shown by the curve made up of dots and labeled al . Simulations were 
carried out with the EPRI-AMPS simulator using a discrete energy state located at the center of the band gap with 
densities of 7E17 and 8E17 cm^eV" 1 . The simulated results are shown by the solid curves in Figure 6 that are 
labeled 7E17 and 8E17. Our previous work shows that there are three current mechanisms that contribute to the 
dark J-V characteristic of irradiated a-Si:H solar cells, namely, injection, shunt and CV m mechanisms [9]. The 
current injection mechanism is evident in the rl curve since the current increases exponentially with forward-bias 
voltage above 0.7 Volts. Below 0.7 V the shunt and CV m mechanisms appear to contribute to the current. The rl 
dark J-V characteristic has higher J 0 and n values than the al characteristic. The simulated results suggest that a 
good fit to the irradiated cell J-V characteristic above 0.7 V may be obtained with a discrete energy state density 
of about 7.5E17 cm^eV' 1 . 

The results of the EPRI-AMPS dark J-V simulations in Figure 5 show that the cell current density for a given 
forward-bias voltage increases with increasing discrete energy state density up to 1 E 1 7 cm^eV' 1 . However, 

Figure 6 shows that increasing the discrete energy state density further to 7E17 and 8E17 cm^eV' 1 results in a 
decrease in the current density at a given voltage. Note that the measured J-V characteristic above 0.7 V shown 
in Figure 6 for the irradiated cell has a lower current density than that of the annealed cell for the same forward- 
bias voltages. Hence, in the forward-bias voltage range where current the injection mechanism dominates, the 
effect of irradiation on the cell is to produce a behavior in the measured J-V characteristic that is the same as 
increasing the discrete energy state density in the simulated dark J-V characteristics. It is plausible to assume 
that 40 keV electron irradiation produces defect states near the center of the band gap that serve as 
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Figure 6. Measured and EPRI-AMPS simulated dark J-V characteristics for an a-Si:H 
p-i-n solar cell before and after irradiation with 40 keV electrons. 


recombination centers. The behavior can also be explained by noting the effect of increasing the discrete energy 
state density, whether by simulation or irradiation, is to increase the values of J 0 and n in Equation (5). The 
simulated dark J-V characteristics plotted in Figure 6 show that the higher discrete energy state densities result in 
a non-linear current density on the log scale below 0.7 V, The same effect exists in the measured dark J-V 
characteristic of the irradiated cell, but it is much more pronounced; the magnitude of the current density is about 
an order of magnitude larger for the measured cells. In our earlier work that employed parametric fitting, we 
associated this effect with a shunting mechanism and a current mechanism that results in a term of the form CV m . 
The EPRI-AMPS simulation studies provide an understanding of the dark J-V characteristic in the forward-bias 
region where the current injection mechanism dominates. Studies need to be carried out to understand the 
physical bases of the current mechanisms at forward biases below 0.7 V. 

CONCLUSIONS 

Studies have been carried out to correlate the effect of 40 keV electron irradiation on dark J-V characteristics. 
EPRI-AMPS and PC-1D simulators were utilized to characterize a two-sided step p-n junction before selecting a 
simulator for investigating a p-i-n solar cell structure. The two simulators produced dark J-V characteristics that 
were in agreement each other and an analytical model, but the lifetimes used in the two simulators differ by a 
factor of ten. The EPRI-AMPS simulator was used to simulate dark J-V characteristics of a p-i-n solar cell. In the 
forward-bias voltage range above 0.7 V, where the current injection mechanism dominates, the effect of 
irradiation on the cell is to produce a behavior in the measured J-V characteristic that is the same as increasing 
the discrete energy state density in the simulated dark J-V characteristics. It is plausible to assume that 40 keV 
electron irradiation produces defect states near the center of the band gap and that the states serve as 
recombination centers. The EPRI-AMPS simulation studies provide an understanding of the dark J-V 
characteristic when current injection dominates. Studies need to be carried out to understand the physical bases 
of the current mechanisms at forward biases below 0.7 V, 
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